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Description 

SYSTEMS, METHODS, AND AN ARTICLE 
OF MANUFACTURE FOR DETERMINING 
FREQUENCY VALUES ASSOCIATED WITH 
FORCES APPLIED TO A DEVICE 

Background of Invention 

[0001] Mechanical systems or devices can be modeled mathe- 
matically using dynamic models. The dynamic models 
generally utilize a forcing waveform as an input for the 
model. The forcing waveform is transferred mathemati- 
cally using a Fourier transform from the time domain to 
the frequency domain. In the frequency domain, the forc- 
ing waveform is represented by a frequency domain spec- 
trum corresponding to a plurality of spectrum lines each 
having a particular amplitude and frequency. Generally, 
the frequency domain spectrum is utilized as an input to 
the dynamic model to compute a desired output wave- 
form. Computing the desired output waveform is also 



called solving the dynamic model. The computational time 
in a computer, however, is relatively high when computing 
the desired output waveform using an entire frequency 
domain spectrum of the forcing waveform. 
[0002] Accordingly, it would be desirable to have a method for 
selecting a subset of the frequency values of the fre- 
quency domain spectrum associated with a forcing wave- 
form in order to reduce the amount of computational time 

required to solve a dynamic model. 
Summary of Invention 

[0003] a method for determining frequency values associated 
with forces applied to a device in accordance with an ex- 
emplary embodiment is provided. The method includes 
determining a first plurality of spectral amplitude values 
associated with a first forcing waveform applied to the 
device. The method further includes determining a second 
plurality of spectral amplitude values associated with a 
second forcing waveform applied to the device. The 
method further includes determining a maximum spectral 
amplitude value based on the first and second plurality of 
spectral amplitude values. The method further includes 
determining a threshold amplitude value based on the 
maximum spectral amplitude value and an acceptance 



value. The method further includes determining a first 
plurality of desired frequency values by selecting fre- 
quency values associated with a subset of the first plural- 
ity of spectral amplitude values that are greater than or 
equal to the threshold amplitude value. Finally, the 
method includes determining a second plurality of desired 
frequency values by selecting frequency values associated 
with a subset of the second plurality of spectral amplitude 
values that are greater than or equal to the threshold am- 
plitude value. 

[0004] An article of manufacture in accordance with an exem- 
plary embodiment is provided. The article of manufacture 
includes a computer storage medium having a computer 
program encoded therein for determining frequency val- 
ues associated with forces applied to a device. The com- 
puter storage medium includes code for determining a 
first plurality of spectral amplitude values associated with 
a first forcing waveform applied to the device. The com- 
puter storage medium further includes code for determin- 
ing a second plurality of spectral amplitude values associ- 
ated with a second forcing waveform applied to the de- 
vice. The computer storage medium further includes code 
for determining a maximum spectral amplitude value 



based on the first and second plurality of spectral ampli- 
tude values. The computer storage medium further in- 
cludes code for determining a threshold amplitude value 
based on the maximum spectral amplitude value and an 
acceptance value. The computer storage medium further 
includes code for determining a first plurality of desired 
frequency values by selecting frequency values associated 
with a subset of the first plurality of spectral amplitude 
values that are greater than or equal to the threshold am- 
plitude value. Finally, the computer storage medium in- 
cludes code for determining a second plurality of desired 
frequency values by selecting frequency values associated 
with a subset of the second plurality of spectral amplitude 
values that are greater than or equal to the threshold am- 
plitude value. 

[0005] a system for determining frequency values associated 

with forces applied to a device in accordance with another 
exemplary embodiment is provided. The system includes 
a first sensor operably coupled to the device. The first 
sensor generates a first signal over time indicative of a 
first forcing waveform applied to the device. The system 
further includes a second sensor operably coupled to the 
device. The second sensor generates a second signal over 



time indicative of a second forcing waveform applied to 
the device. The system further includes a computer oper- 
ably communicating with the first and second sensors. 
The computer is configured to determine a first plurality 
of spectral amplitude values associated with the first forc- 
ing waveform. The computer is further configured to de- 
termine a second plurality of spectral amplitude values 
associated with the second forcing waveform. The com- 
puter is further configured to determine a maximum 
spectral amplitude value based on the first and second 
plurality of spectral amplitude values. The computer is 
further configured to determine a threshold amplitude 
value based on the maximum spectral amplitude value 
and an acceptance value. The computer is further config- 
ured to determine a first plurality of desired frequency 
values by selecting frequency values associated with a 
subset of the first plurality of spectral amplitude values 
that are greater than or equal to the threshold amplitude 
value. The computer is further configured to determine a 
second plurality of desired frequency values by selecting 
frequency values associated with a subset of the second 
plurality of spectral amplitude values that are greater than 
or equal to the threshold amplitude value. 



[0006] a system for determining frequency values associated 

with forces applied to a device in accordance with another 
exemplary embodiment is provided. The system includes 
a first sensor means operably coupled to the device for 
generating a first signal over time indicative of a first 
forcing waveform applied to the device. The system fur- 
ther includes a second sensor means operably coupled to 
the device for generating a second signal over time in- 
dicative of a second forcing waveform applied to the de- 
vice. The system further includes a computer means for 
operably communicating with the first and second sen- 
sors. The computer means is configured to determine a 
first plurality of spectral amplitude values associated with 
the first forcing waveform. The computer means is further 
configured to determine a second plurality of spectral am- 
plitude values associated with the second forcing wave- 
form. The computer means is further configured to deter- 
mine a maximum spectral amplitude value based on the 
first and second plurality of spectral amplitude values. 
The computer means is further configured to determine a 
threshold amplitude value based on the maximum spec- 
tral amplitude value and an acceptance value. The com- 
puter means is further configured to determine a first 



plurality of desired frequency values by selecting fre- 
quency values associated with a subset of the first plural- 
ity of spectral amplitude values that are greater than or 
equal to the threshold amplitude value. The computer 
means is further configured to determine a second plural- 
ity of desired frequency values by selecting frequency val- 
ues associated with a subset of the second plurality of 
spectral amplitude values that are greater than or equal to 
the threshold amplitude value. 
[0007] other systems and/or methods according to the embodi- 
ments will become or are apparent to one with skill in the 
art upon review of the following drawings and detailed 
description. It is intended that all such additional systems 
and methods be within the scope of the present invention, 

and be protected by the accompanying claims. 
Brief Description of Drawings 

[0008] Figure 1 illustrates a compressor monitoring system in- 
cluding a compressor and a monitoring computer. 

[0009] Figure 2 is a schematic of a portion of a crankshaft of the 
compressor of Figure 1. 

[0010] Figures 3-8 are flowcharts of a method for determining 
frequency values associated with forces applied to the 
compressor crankshaft. 



[0011] Figure 9 is a simplified schematic of the crankshaft, the 

connecting rod, the piston rod, and the piston — of the 

compressor illustrated in Figure 1. 
Detailed Description 

[0012] Referring to Figures 1, 2 and 9, a compressor monitoring 
system 10 is illustrated for monitoring the operation of a 
compressor 12. The compressor monitoring system 10 
further includes pressure sensors 15, 16, 17, 18 and a 
monitoring computer 14. The configuration of the com- 
pressor monitoring system 10 will be briefly explained in 
order to better understand a method for determining fre- 
quency values associated with forces applied to a 
crankshaft of the compressor 12. 

[0013] The compressor 12 is provided to force a fluid through a 
head-end chamber 32 and a crank end chamber 34. As 
shown, compressor 12 includes a housing 19, a 
crankshaft 20, a connecting rod 25, a cross-head 26, a 
cross-head pin 124, a piston rod 28, and a piston 30. It 
should be noted that although only one piston is shown 
for simplicity of understanding, compressor 12 can in- 
clude a plurality of pistons, connecting rods, cross-head 
pins, and piston rods. 

[0014] Th e housing is provided to enclose all of the remaining 



elements of the compressor 12. The crankshaft 20 is dis- 
posed in the housing 19 and is operably coupled to an ac- 
tuating means (not shown) that rotates crankshaft 20. The 
crankshaft 20 has a first crankshaft throw 21, adjacent a 
large end 120 of connecting rod 25. Further, the 
crankshaft 20 is operably coupled to the large end 120 of 
connecting rod 25. Further, the cross-head 26 is operably 
coupled at the cross-head pin 124 to a small end 122 of 
connecting rod 25. 

[0015] it should be noted that crankshaft 20 further includes a 
crankshaft throw 22 torsionally coupled adjacent another 
connecting rod (not shown). 

[0016] As shown, a second end of the cross-head 26 is coupled 
to the piston rod 28 which is linearly driven along an axis 
35 by the cross-head 26. The piston rod 28 is further 
coupled to the piston 30. Thus, linear movement of the 
piston rod 28 linearly moves the piston 30 into either 
head-end chamber 32 or crank-end chamber 34, depend- 
ing upon the direction of movement of the piston 30. 

[0017] The pressure sensors 15, 16 are provided to generate 
pressure signals (PI) and (P2) respectively, indicative of 
the pressure generated by piston 30 in the crank-end 
chamber 34 and the head-end chamber 32 respectively. 



As will be explained in greater detail below, signals (PI) 
and (P2) will be utilized by computer 14 to determine a 
forcing waveform representing forces applied to a plane 
24 of the crankshaft throw portion 21. 

[0018] The pressure sensors 17, 18 are provided to generate 
pressure signals (P3) and (P4) respectively, indicative of 
the pressure generated by another piston (not shown) in a 
crank-end chamber and a head-end chamber. As will be 
explained in greater detail below, signals (P3) and (P4) will 
be utilized by computer 14 to determine a forcing wave- 
form representing forces applied to a plane 23 of the 
crankshaft throw portion 22. 

[0019] The monitoring computer 14 is provided to receive the 

pressure signals (PI), (P2), from pressure sensors 15, 16, 
respectively and to generate a first forcing waveform. The 
computer 14 is further provided to receive the pressure 
signals (P3), (P4) from the pressure sensors 17, 18 and to 
generate a second forcing waveform. After generating the 
first and second forcing waveforms, the computer 14 is 
provided to implement a method for determining fre- 
quency values associated with the first and second forcing 
waveforms 20. The monitoring computer 14 includes a 
CPU 36 that operably communicates with the storage me- 



dia including read-only memory (ROM) 38 and a random 
access memory (RAM) 40. The storage media may be im- 
plemented using any of a number of known memory de- 
vices such as PROMs, EPROMs, EEPROMS, flash memory or 
any other electric, magnetic, optical or combination mem- 
ory device capable of storing data, some of which repre- 
sent executable instructions used by CPU 36. The CPU 36 
communicates via the I/O interface 42 with pressure sen- 
sors 15, 16, 17, 18. 
[0020] Referring to Figures 3-8, a method for determining fre- 
quency values associated with forces applied to the 
crankshaft 20 will now be explained. The method is di- 
rected to determining a first and second plurality of de- 
sired frequency values associated with first and second 
forcing waveforms, respectively, applied to the crankshaft 
throw portions 21, 22, respectively. It should be noted 
that additional forcing waveforms could be calculated for 
additional components associated with the crankshaft 20 
including crankshaft throw portions, a flywheel, and a 
motor driving the crankshaft (not shown). Further, it 
should be noted that in the exemplary embodiment the 
CPU 36 can execute the steps 60-78 utilizing a first forc- 
ing waveform while concurrently executing steps 80-98 



utilizing a second forcing waveform. 
[0021] At step 60, the CPU 36 determines a first gas force wave- 
form relating to a gas force acting on the first piston 30 
based on the pressure signals (PI), (P2). The first gas 
force waveform can be determined by iteratively calculat- 
ing the following equation over time: 

Fgas) = {PcfaAcex — PhexAhex) — PambiArodx 
[0022] where 

[0023] f represents an instantaneous gas pressure; 

GASl 

[0024] p represents the pressure in the crank-end chamber 

obtained from the signal PI; 
[0025] a cei represents the area of the crank-end chamber; 

[0026] p represents the pressure in the head-end chamber 

obtained from the signal P2; 
[0027] a represents the area of the head-end chamber 

HEl r 

[0028] P ambi represents the pressure of gas in the housing push- 
ing against the piston rod; and 
[0029] a represents the area of the piston rod. 

RODl 

[0030] At step 62, the CPU 36 determines a first inertia force 

waveform relating to an inertia force of the first cylinder 



piston 30. The first inertia force waveform can be deter- 
mined by iteratively calculating the following equation 
over time: 

Fm\ — (fftxheadl + Tttpl + Hip rod] + TYlconsmX )xl 
[0031] W here 

[0032] p represents an instantaneous inertia force of the cross- 
mi 

head pin, the piston, the piston rod, and the connecting 
rod small end; 

[0033] m represents the mass of the cross-head pin; 

xheadl K K 

[0034] m ^ represents the mass of the piston; 

[0035] m^ odi represents the mass of the piston rod; and, 

[0036] m represents the mass of the connecting rod small 

consml 

end 

[0037] ■ - 

X 

1 represents the acceleration of the cumulative masses il- 
lustrated in the equation. 
[0038] At step 64, the CPU 36 determines a first forcing wave- 
form indicative of a crankshaft torque at the crankshaft 
throw 21. The first forcing waveform is determined based 



on the first gas force waveform and the first inertia force 
waveform. The first forcing waveform can be determined 
by iteratively calculating the following equation over time: 

71 = -{Fgasi + Fm\)r sin/91 
[0039] where 

[0040] ti represents an instantaneous torque at the crankshaft 
throw; 

[0041] r represents the distance from a crankshaft rotation axis 

127 to the crankshaft yoke centerpoint 126; 
[0042] 0i represents an angular position of the crankshaft; and, 

[0043] |_ represents the length of the connecting rod between the 
crankshaft yoke centerpoint 126 and the centerpoint of 
the cross-head pin 124. 

[0044] N ext at step t he CPU 36 makes a determination as to 
whether the first forcing waveform was generated over an 
integral number of revolutions of the crankshaft. If the 
value of step 66 equals "yes", the method advances to 
step 70. Otherwise the method advances to step 68. 

[0045] At step 68, the CPU 36 copies portions of the first forcing 
waveform to itself to obtain a first forcing waveform hav- 



1 + 



r cos 6\ 



Vi 2 - r 1 sin 2 01 



ing points relating to an integral number of revolutions of 
the crankshaft. After step 68, the method advances to 
step 70. 

[0046] At step 70, the CPU 36 removes a first DC component 

from the first forcing waveform to obtain a first modified 
forcing waveform. 

[0047] At step 72, CPU 36 stores the first DC component in ROM 
38. 

[0048] At step 74, the CPU 36 applies a Fourier transform to the 
first modified forcing waveform to obtain a first plurality 
of complex spectral values. For example, the first plurality 
of complex spectral values could have the following val- 
ues: (i) 1.2+jl.6, (ii) 1.0+jl.5, (iii) 1.5+j2.1, (iv) 0.8+j0.2. 

[0049] N ext at step jq^ t he CPU 36 calculates a first plurality of 
spectral amplitude values based on the first plurality of 
complex spectral values. Each of the first plurality of 
spectral amplitude values can be calculated using the fol- 
lowing equation: 

Ampl = VRe 2 + Im 2 

[0050] W here 



[0051] Ampl represents the spectral amplitude; 



[0052] 




represents the square of the real portion of the complex 
number; 



[0053] 




represents the square of the imaginary portion of the 
complex number. 
[0054] For example, the first plurality spectral amplitude values 
could have the values: (i) 2.0, (ii) 1.8, (iii) 2.6, (iv) 0.8 — 
based on the first plurality of complex spectral values of: 
(i) 1.2+11.6, (ii) 1.0+jl.5, (iii) 1.5+j2.1, (iv) 0.8+j0.2, re- 
spectively. 

[0055] At step 78, the CPU 36 determines a first maximum spec- 
tral amplitude from the first plurality of spectral amplitude 
values. In particular, the CPU 36 determines which one of 
the first plurality of spectral amplitude values has the 
greatest numerical value. 

[0056] Referring to Figures 1, 3, 6, and 7, the steps 80-98 gen- 
erate a second forcing waveform associated with a second 
piston (not shown) that will now be explained. For pur- 
poses of discussion, the second piston has an associated 



crank-end chamber, head-end chamber, piston, cross- 
head pin, and connecting rod. 
[0057] At step 80, the CPU 36 determines a second gas force 

waveform relating to a gas force acting on a second piston 
(not shown) based on the pressure signals (P3), (P4). It 
should be noted that the sampling start time of pressure 
signals (P3) and (P4) is identical to the sampling start time 
of pressure signals (PI) and (P2). Further, the sampling 
rate of pressure signals (P3) and (P4) is identical to the 
sampling rate of pressure signals (PI) and (P2). The sec- 
ond gas force waveform can be determined by iteratively 
calculating the following equation over time: 

FOAS2 = (PcE2AcE2-PhE2AhE2)-PaMB2ArOD2 



[0058] W here 

[0059] p represents an instantaneous gas pressure; 

GAS 2 

[0060] p represents the pressure in the crank-end chamber 

obtained from the signal P3; 
[0061] a represents the area of the crank-end chamber; 

[0062] p represents the pressure in the head-end chamber 

obtained from the signal P4; 
[0063] a represents the area of the head-end chamber; 

HE2 r 



[0064] P amb2 represents the pressure of gas in the housing push- 
ing against the piston rod; and 
[0065] a represents the area of the piston rod. 

ROD2 K K 

[0066] At step 82, the CPU 36 determines a second inertia force 
waveform relating to an inertia force of the second cylin- 
der piston (not shown). The second inertia force waveform 
can be determined by iteratively calculating the following 
equation over time: 

Fm2 = (nixheadl + Mpl + Jflprodl + M consm 

[0067] W here 

[0068] p represents an instantaneous inertia force of the cross- 
head pin, the piston, the piston rod, and the connecting 
rod small end; 

[0069] m represents the mass of the cross-head pin; 

xhead2 K ^ 

[0070] m represents the mass of the piston; 

[0071] m represents the mass of the piston rod; and, 

prod 2 

[0072] m represents the mass of the connecting rod small 

consm2 



[0073] 



end 



x 



2 represents the acceleration of the cumulative masses il- 
lustrated in the forgoing equation. 
[0074] At step 84, the CPU 36 determines a second forcing wave- 
form indicative of a crankshaft torque at the crankshaft 
throw 22. The second forcing waveform is determined 
based on the second gas force waveform and the second 
inertia force waveform. The second forcing waveform can 
be determined by iteratively calculating the following 
equation over time: 

rcos#2 
4l 2 -r 2 sin 2 6>2„ 

[0075] W here 

[0076] j2 represents an instantaneous torque at the crankshaft 
throw; 

[0077] r represents the distance from a crankshaft rotation axis 

127 to the crankshaft yoke centerpoint; 
[0078] 02 represents an angular position of the crankshaft; 

[0079] |_ represents the length of the connecting rod between the 
crankshaft yoke centerpoint and the centerpoint of the 
cross-head pin. 

[0080] N ex t at step 86, the CPU 36 makes a determination as to 



T2 = -{Fgasi + Fmi)r sin 02 



whether the second forcing waveform was generated over 
an integral number of revolutions of the crankshaft. If the 
value of step 86 equals "yes", the method advances to 
step 90. Otherwise, the method advances to step 88. 

[0081] At step 88, the CPU 36 copies portions of the second 
forcing waveform to itself to obtain a second forcing 
waveform having points relating to an integral number of 
revolutions of the crankshaft. After step 88, the method 
advances to step 90. 

[0082] At step 90, the CPU 36 removes a second DC component 
from the second forcing waveform to obtain a second 
modified forcing waveform. 

[0083] At step 92, CPU 36 stores the second DC component in 
ROM 38. 

[0084] At step 94, the CPU 36 applies a Fourier transform to the 
second modified forcing waveform to obtain a second 
plurality of complex spectral values. For example, the 
second plurality of complex spectral values could have the 
following values: (i) 1.2+j 1.6, (ii) 1.0+jl.5, (iii) 1.5+j2.1, 
(iv) 0.8+j0.2. 

[0085] Next to step 96, the CPU 36 calculates a second plurality 
of spectral amplitude values based on the second plurality 
of complex spectral values. Each of the second plurality of 



spectral amplitude values can be calculated using the fol- 
lowing equation: 




[0086] W here 



[0087] Amp2 represents a spectral amplitude; 



[0088] 



Re 2 



represents the square of the real portion of the complex 
number; 



represents the square of the imaginary portion of the 
complex number. 

[0090] For example, the second plurality spectral amplitude val- 
ues could have the values: (i) 2.0, (ii) 1.8, (iii) 2.6, (iv) 0.8 
— based on the second plurality of complex spectral val- 
ues of: (i) 1.2+J1.6, (ii) 1.0+jl.5, (iii) 1.5+j2.1, (iv) 
0.8+j0.2, respectively. 

[0091] At step 98, the CPU 36 determines a second maximum 



[0089] 




spectral amplitude from the second plurality of spectral 
amplitude values. In particular, the CPU 36 determines 
which one of the second plurality of spectral amplitude 
values has the greatest numerical value. 
[0092] Referring to Figures 5 and 8, after either of steps 78 or 
98, the method advances to step 100. At step 100, the 
CPU 36 determines an overall maximum spectral ampli- 
tude by calculating the greater of the first maximum 
spectral amplitude and the second maximum spectral am- 
plitude. 

[0093] Next to step 102, the CPU 36 determines a threshold am- 
plitude value based on the overall maximum spectral am- 
plitude and acceptance of value. For example, the thresh- 
old amplitude value can be calculated using the following 
equation: 

[0094] threshold amplitude value = overall maximum spectral 
amplitude * acceptance value 

[0095] The acceptance value can be empirically determined based 
upon the number of desired spectral amplitudes and the 
desired degree of accuracy in the solution of the dynamic 
model. For example, the overall maximum spectral ampli- 
tude could be 2.6 and the acceptance value could be 
about 0.4, resulting in a threshold amplitude value of 



1.04. It should be noted that the desired degree of accu- 
racy can be empirically determined by comparing the so- 
lution of the dynamic model using a specific acceptance 
value to a corresponding measured value. 

[0096] Next at step 104, the CPU 36 determines a first plurality 
of desired frequency values by selecting frequency values 
associated with a subset of the first plurality of spectral 
amplitude values that are greater than or equal to the 
threshold amplitude value. For example, the first plurality 
of desired frequency values could include: frequencyl, 
frequency2, and frequency3 because the spectral ampli- 
tude values associated with these frequencies are greater 
than the threshold amplitude value of 1.04. 

[0097] Next at step 106, the CPU 36 determines a second plural- 
ity of desired frequency values by selecting frequency val- 
ues associated with a subset of the second plurality of 
spectral amplitude values that are greater than or equal to 
the threshold amplitude value. For example, the second 
plurality of desired frequency values could include: fre- 
quency5, frequency6, and frequency7 because the spec- 
tral amplitude values associated with these frequencies 
are greater than the threshold amplitude value of 1.04. 

[0098] Next at step 108, the CPU 36 solves a dynamic model of 



the crankshaft using a set of frequencies that comprise 
the union of: (i) the first and second DC components, (ii) a 
subset of the first plurality of complex spectral values as- 
sociated with the first plurality of desired frequency val- 
ues, and (iii) a subset of the second plurality of complex 
spectral values associated with the second plurality of de- 
sired frequency values. 

[0099] The inventive method and article of manufacture for de- 
termining frequency values associated with forces applied 
to the crankshaft represents a substantial advantage over 
known methods. In particular, the embodiments of the in- 
vention provide a technical effect of selecting a subset of 
the frequency values of the frequency domain spectrum 
associated with one or more forcing waveforms in order to 
dramatically reduce the amount of computational time re- 
quired to solve a dynamic model in a computer. 

[0100] while the invention is described with reference to an ex- 
emplary embodiment, it will be understood by those 
skilled in the art that various changes may be made and 
equivalence may be substituted for elements thereof with- 
out departing from the scope of the invention. In addition, 
many modifications may be made to the teachings of the 
invention to adapt to a particular situation without de- 



parting from the scope thereof. Therefore, it is intended 
that the invention not be limited to the embodiment dis- 
closed for carrying out this invention, but that the inven- 
tion includes all embodiments falling with the scope of 
the intended claims. Moreover, the use of the term's first, 
second, etc. does not denote any order of importance, but 
rather the term's first, second, etc. are used to distinguish 
one element from another. 



